The influence of urban environmental effects on the orchard soil microbial community structure and function: a case study in Zhejiang, China

The urban environmental effects can have multifaceted impacts on the orchard soil microbial community structure and function. To specifically study these effects, we investigated the soil bacterial and fungal community in the laxly managed citrus orchards using amplicon sequencing. Ascomycota demonstrated significant dominance within the citrus orchard soils. The increased presence of beneficial Trichoderma spp. (0.3%) could help suppress plant pathogens, while the elevated abundance of potential pathogenic fungi, such as Fusarium spp. (0.4%), might raise the likelihood of disorders like root rot, thereby hindering plant growth and resulting in reduced yield. Moreover, we observed significant differences in the alpha and beta diversity of bacterial communities between urban and rural soils (p < 0.001). Environmental surveys and functional prediction of bacterial communities suggested that urban transportation factors and rural waste pollution were likely contributing to these disparities. When comparing bacterial species in urban and rural soils, Bacillus spp. exhibited notable increases in urban areas. Bacillus spp. possess heavy metal tolerance attributed to the presence of chromium reductase and nitroreductase enzymes involved in the chromium (VI) reduction pathway. Our findings have shed light on the intricate interplay of urban environmental effects and root systems, both of which exert influence on the soil microbiota. Apart from the removal of specific pollutants, the application of Bacillus spp. to alleviate traffic pollution, and the use of Trichoderma spp. for plant pathogen suppression were considered viable solutions. The knowledge acquired from this study can be employed to optimize agricultural practices, augment citrus productivity, and foster sustainable agriculture.


Introduction
Due to global climate change, microorganisms impacting human health and societal productivity underwent changes.This included the emergence of new pathogens (Casadevall, 2023), alterations in their growth and reproductive characteristics, particularly the increase in antibiotic resistance (Lockhart et al., 2023), the expansion of their distribution range, and seasonal fluctuations in microbial infectious diseases (Mora et al., 2022).The urban environment served as the primary habitat for modern human populations and played a crucial role in the interaction between humans and nature.For local regions, urban environmental effects such as the urban heat island effect and pollution emissions could have even more direct impacts on microorganisms than the effects of global climate warming (Parajuli et al., 2018;Abrego et al., 2020;Li et al., 2023).Toxic heavy metals, organic pollutants, emerging contaminants, and other biotic and abiotic stressors may impact nutrient utilization, plant metabolic pathways, agricultural productivity, and soil fertility (Pathak et al., 2024).Therefore, investigating the effects of urban environmental effects on microbial community structure could provide valuable insights into the distribution and transmission patterns of microorganisms and facilitate the assessment of environmental health risks (Blocker et al., 2020;Zhou and Zhou, 2023).
The soil environment serves as a reservoir for pathogenic microorganisms and is also a critical medium linking microbial communities and human activities (Williams et al., 2024;Wang B. et al., 2024;Wang X. et al., 2024).Pathogenic microorganisms in the soil not only directly impacted agricultural production but also resulted in alterations in the human microenvironment and even clinical infections (Liu et al., 2024;Singh et al., 2024;Yiallouris et al., 2024).In fact, soil microorganisms play crucial roles in shaping soil health (Cheng et al., 2019;Li et al., 2022), facilitating nutrient cycling (Hu et al., 2021), and fostering the overall health of the ecosystem (Tsitsigiannis et al., 2005).They participate in various biogeochemical processes, such as decomposition, nitrogen fixation (Pankievicz et al., 2021;Yin et al., 2021), and organic matter recycling (Yang et al., 2022).Moreover, microbes in rhizospheric soil engage in symbiotic interactions with plant roots (Matilla et al., 2007), influencing nutrient uptake (Wu et al., 2022), stress tolerance (Igiehon et al., 2021), and disease resistance (De Tender et al., 2021).Understanding the structure and function of microbial communities, particularly in the rhizospheric soil, is crucial for comprehending the distribution and transmission patterns of microorganisms, promoting agricultural production, and assessing the environmental health risks posed by urban effects (Aguilera et al., 2022;Furlan et al., 2023;Pastrana et al., 2023).Research on citrus has shown that soil microbial communities, especially those associated with the root system, significantly influence the quality of citrus fruits and interact with the host immune system (Su et al., 2023).
The warm and humid climate of Zhejiang Province, China, characterized by hot summers and mild winters, provides an ideal environment for the cultivation of citrus fruits, establishing citrus as a traditional fruit in the region (Lin et al., 2023).Given the lax soil management practices observed in this area, where active managerial intervention is infrequent, there emerges an opportunity to delve into the effects of urban environmental factors on soil microbiota within authentic field conditions.This was because the Zhejiang Province had a developed urban economy, with the primary, secondary, and tertiary industries accounting for 2.9, 41.9, and 55. 2%, respectively (Jinhua City Statistics Bureau, 2023).Labor in rural areas tended to concentrate toward urban centers, leading plenty of older farmers to preferentially choose easily managed citrus for cultivation.Therefore, this study took citrus orchard soil environments close to their native state as examples.Soil samples were collected from the rhizospheric and peripheral soil at different distances from urban areas during the fruit ripening period on family farms.The microbial community structure and functional profiles in the soil were analyzed using sequencing methods targeting the 16S rRNA gene and ITS region gene.
The fruit ripening period is one of the most frequent times for farmers to interact with the soil and marks a critical phase in the life cycle of citrus trees, characterized by intricate physiological and ecological changes.Metabolism during this period is a key determinant of citrus flavor and nutritional quality (Saini et al., 2020;Sheng et al., 2022), and troubles during fruit ripening can lead to significant production losses (Olimi et al., 2022;Sharma et al., 2023;Kifle et al., 2024).Additionally, during this phase, the citrus orchard ecosystem experiences more frequent interactions, including proliferation of saprophytic microorganisms due to fruit rot and the presence of excreta from animals, primarily birds (Thompson and Willson, 1979;Schaefer and Ruxton, 2011;Whitehead and Poveda, 2011).Therefore, during this period, the impacts of urban environmental effects such as noise pollution and vehicular movement are more sensitively captured (Moore et al., 2002;Lopez-Velasco et al., 2012;Xi et al., 2016).
The objective of this study is to assess the impact of urban environmental effects on microbial distribution and environmental health risks through the compositional differences of microbes in the soil of citrus orchards in their native state during the fruit ripening period.Furthermore, the study aims to provide valuable evidence for research on environmental-host relationships, optimization of future agricultural production, and monitoring of the prevalence of pathogenic microorganisms using these research findings and microbial community data.

Study site
The soil sampling was carried out in Jinhua, situated in the central region of Zhejiang Province, China.The sampling commenced on October 9, 2022, and continued for 10 days, with each sampling site's process being completed within a single day.Jinhua is located approximately between 119°14′-120°47′ east longitude and 28°32′-29°41′ north latitude, spanning 129 km north to south and 151 km east to west.The orchard area covers 26,290 hectares, with 3,809 hectares dedicated to citrus orchards.This region falls within the subtropical zone, characterized by a warm and humid climate, with hot summers and mild winters.In 2022, the total annual sunshine hours were 1701.2, with a precipitation of 1300.3 mm and an average temperature of 19.4°C.
To ensure the similarity of sampling environments and eliminate other factors' influence on soil microorganisms, we specifically selected family farms adjacent to rivers for environmental sampling.Four villages with lax management were chosen as sampling sites based on their geographic location.Lax management is defined as a management practice where, apart from essential watering, the intervention frequency (including pesticide and fertilizer usage) is equal to or less than twice a year; family farms with higher intervention frequencies were not considered.Among these, priority was given to farms that either abandoned management or had lower management frequencies.
Village A located on the river downstream of the center of the city and closely adjacent to the city.It is situated in a densely populated area and is strongly influenced by urban population mobility.On the other hand, villages B, C, and D are situated farther from the city and are next to the river but upstream of the city (Figure 1A).Fourteen family farms were designated as individual sampling site (Supplementary Table S1), where five sampling points were selected: one sampling point was in the center of the farm, and one was at each corner.In order to assess disparities in composition of the soil in the rhizosphere and the soil outside of the rhizosphere, two additional sampling points were strategically placed outside the rhizosphere (peripheral soil).Specifically, these points were positioned at the midpoint between two fruit trees, where no root growth was observed.

Sample collection
We conducted sampling using five-point sampling method during the ripening stages of Jinhua's native citrus.Rhizospheric soil samples  were gathered from the rhizosphere, which was considered to be the area around the root system.Five citrus trees aged 8-10 years were selected at each sampling point, and soil samples were collected around the roots of each tree to a depth of approximately 20 cm.We did not sample soil from the roots of citrus trees showing signs of disease or root damage.For the two sampling points located at the midpoint between two fruit trees, peripheral soil samples were collected to the same depth (20 cm), but without including the rhizospheric soil.During soil excavation, areas within a 1-m radius of the soil surface containing evident contaminants and decayed matter were excluded.
Before sampling, we removed the topsoil using a shovel.For sampling of rhizospheric soil, we gently removed soil that was loosely attached to the roots.Every 10 mL of soil from the five locations within each sampling point was placed in separate sterile plastic bags and thoroughly mixed to create a homogenized composite sample.Approximately 15 mL of the homogenized soil was then transferred into a labeled centrifuge tube and stored at −20°C during transportation.Upon arrival at the laboratory, the samples were promptly transferred to a −80°C ultra-low temperature freezer for long-term storage.The detailed locations of the sampling points and relevant environmental information can be found in Supplementary Table S1.
Finally, we systematically collected a total of 98 soil samples from four villages.Within this dataset, villages A (urban), C (rural), and D (rural) were each represented by 21 samples.Village B (rural) contributed a comprehensive set of 35 soil samples due to its broader and more diverse citrus cultivation landscape.

Extraction of fungal and bacterial DNA
DNA extraction from soil samples was carried out using the QIAamp DNA Mini Kit (Qiagen, Hilden, Germany) following the manufacturer's instructions.To assess the quality and concentration of the extracted DNA, a spectrophotometer (Nanodrop 2000; Thermo Fisher Scientific, Waltham, MA, United States) was employed.The extracted DNA was then stored at −80°C until PCR analysis.
The first round of PCR amplification was performed using the following cycling conditions: 3 min at 95°C; 25 cycles of 30 s at 95°C, 30 s at 55°C, and 30 s at 72°C; and a final extension step at 72°C for 5 min.Each 25 μL PCR mixture contained 12.5 μL of 2× KAPA HiFi HotStart ReadyMix, 1 μL of each primer (1 μM), and 12.5 ng of template DNA.Post-PCR purification was achieved using AMPure XP beads.
The second round of PCR amplification was performed using the following cycling conditions: 3 min at 95°C; 8 cycles of 30 s at 95°C, 30 s at 55°C and 30 s at 72°C; and a final extension at 72°C for 5 min.Each 50 μL PCR mixture contained 25 μL of 2× KAPA HiFi HotStart ReadyMix, 5 μL of each primer, 10 μL of water, and 5 μL of template DNA from products purified in the previous step.A second round of AMPure XP beads purification was performed.
Following purification, the amplicons were pooled in equimolar proportions and subjected to paired-end sequencing on an Illumina MiSeq platform (Illumina, San Diego, United States), according to the standard protocol provided by Majorbio Bio-Pharm Technology Co. Ltd.The number of sequencing reads obtained can be found in Supplementary Table S2.The sequencing reads have been deposited in the NCBI BioProject under ID PRJNA1007597.
Taxonomy assignment was performed on all representative sequences after training the species classifier using the q2-featureclassifier (Bokulich et al., 2018) plugin.The training process involved utilizing primers to extract target sequences, followed by filtering out corresponding taxonomic information using the RESCRIPt (Robeson et al., 2021) plugin.Finally, the feature-classifier plugin was used for training.The results of taxonomic analyses indicated that 98.0% of the bacteria were successfully identified at the phylum level using the 16S rRNA gene, whereas only 67.6% of fungi were successfully identified.Visualization of the taxonomy results was carried out on QIIME2 View 1 .The amplicon package 2 of R was used to illustrate species compositions and to create Venn diagrams and circle plots.
Functional prediction of the 16S rRNA gene representative sequences was performed with the picrust2 (Douglas et al., 2020) plugin within the QIIME2 pipeline.For BugBase (Ward et al., 2017) phenotypic predictions, submissions were made via the online platform 3 .The Linear Discriminant Analysis Effect Size (LEfSe) algorithm was processed using the format2lefse function within the amplicon package of R, followed by submission to ImageGP (Keighley et al., 2022)  For functional classification of the ITS region gene, FUNGuild (Nguyen et al., 2016) was employed to parse fungal community datasets based on trophic mode, trait and growth form.Differences between two groups were performed with the Mann-Whitney U test.

Diversity and statistical analysis
Representative sequences were processed using the EasyAmplicon (Liu Y. X. et al., 2021) pipeline.Normalization was conducted using the "otutab_rare.R" script 4 with default parameters.The parameter "--depth 0" was utilized to automatically determine the minimum rarefaction depth.The diversity plugin of QIIME2 was used to calculate alpha and beta diversity.Alpha diversity was assessed using the Abundance-based Coverage Estimator (ACE) index, and differences were compared using the Kruskal-Wallis test.Principal Co-ordinates Analysis (PCoA) based on Bray-Curtis distance was performed.The resulting data were visualized using the amplicon package and pheatmap package of R. 5 Differential analysis was conducted using "EdgeR" and results were visualized by using the 'compare_volcano.R' 6 , 'compare_manhattan.sh'(), and 'compare_ heatmap.sh' 7scripts to generate volcano plots, Manhattan plots, and comparative heatmaps, respectively.A significance threshold of p < 0.05 was considered for all statistical analyses.

Structure and function of microbial communities in urban and rural citrus orchard soils
Among the successfully identified bacteria, the top seven abundant phyla were Proteobacteria, Acidobacteria, Actinobacteria, Chloroflexi, AD3, Gemmatimonadetes, and Firmicutes (Figure 1B).For fungi, the top seven abundant phyla were Ascomycota, Basidiomycota, Chytridiomycota, Mortierellomycota, Rozellomycota, Glomeromycota, and Blastocladiomycota.Notably, Ascomycota fungi were predominant across almost all of the soil samples, accounting for an average of 66.9% of the species (Figure 1C).
Rarefaction curve analyses of the sequencing results demonstrated stable patterns across all samples.From the rarefaction curves, it appeared that village A, which is urban, might have a significantly higher bacterial abundance as compared to the other three villages, which are considered rural (Figure 2A).These intriguing insights were confirmed by comparisons of the alpha diversity values.With regard to pairwise comparisons of bacterial alpha diversity, no statistically significant differences were observed between villages B, C, and D. In contrast, pairwise comparisons of the bacterial alpha diversity values 4 github.com/YongxinLiu/EasyMicrobiome/blob/master/script/otutab_rare.R 5 github.com/cran/pheatmap6 github.com/YongxinLiu/EasyMicrobiome/blob/master/script/compare_volcano.R 7 github.com/YongxinLiu/EasyMicrobiome/blob/master/script/compare_heatmap.shfrom villages A identified significant differences with the alpha diversity values from villages B (p < 0.001), C (p < 0.001), and D (p < 0.001) (Figure 2C and Table 1).We conducted PCoA analysis using pairwise comparisons of Bray-Curtis distances between samples.This analysis confirmed the significant bacterial dissimilarity between urban soil and rural soil (Figure 2E).Furthermore, at the phylum taxonomic level, a noticeable increase in species from Acidobacteria-6 and decreases in Actinobacteria and Ktedonobacteria species were observed in urban soil (Figure 2F).
The most prevalent fungal genera in the samples were Talaromyces, Knufia, Fusarium, Chaetomium, Penicillium, Aspergillus, Paracamarosporium, Trichoderma, and Chaetomella.In terms of fungal alpha diversity, there were no statistically significant differences observed among the villages collectively (p = 0.787) (Figures 2B,D).The results of PCoA analysis revealed only subtle differences between urban soil and rural soil (Supplementary Figure S1A), indicating a generally similar fungal community structure between these two soil environments.

Structure and function of microbial communities in rhizospheric and peripheral citrus orchard soils
In Figures 3A,B, we quantified highly abundant bacteria and fungi (with average relative abundances exceeding 0.2%) using maximum likelihood trees to depict the bacterial and fungal community structures and phylogenetic relationships within the rhizospheric and peripheral soils.This approach provided insights into the community composition and evolutionary connections of bacteria and fungi in these distinct soil environments.
Within the rhizospheric zone, noticeable increases were observed in six bacterial taxa, including Xanthomonadales, Bacilli, and Actinobacteria (Figures 4A,C).Similarly, 27 fungal taxa exhibited heightened representation, encompassing Agaricomycetes, Sordariomycetes, and Eurotiomycetes (Figure 4B).To further elucidate the effects of urban environmental factors and citrus root systems on soil microbiota, we classified samples simultaneously based on different sampling locations and soil types (Supplementary Figure S2).In villages B, C, and D, compared to peripheral soil, the alpha and beta diversity of rhizospheric soil exhibited greater fluctuations, the trend observed in both bacteria and fungi.Specifically, we observed that rhizospheric alpha diversity from the same village had larger interquartile ranges in boxplots compared to peripheral soil.Correspondingly, in the PCA plots of beta diversity, the confidence ellipsoids of rhizospheric samples from the same locations were more inclined to encompass those of peripheral soil.Conversely, in village A, the boxplots of alpha diversity for both bacteria and fungi indicated larger interquartile ranges in peripheral soil.Simultaneously, in the PCA plots, the confidence ellipsoids of peripheral soil and rhizosphere samples exhibited roughly equal intersection.These findings suggest that citrus root systems do indeed play a role in shaping soil microbial community structure, albeit weaker than the influence of urban environmental effects.Furthermore, we contrasted bacterial functional disparities between citrus orchard rhizospheric and peripheral soils (including both urban and rural), pinpointing a significant reduction in the CMP-pseudaminate biosynthesis pathway in the rhizospheric soil (p = 0.001) (Figure 5A).Importantly, the divergence of this pathway was unaffected by urban factors, as evidenced by its substantial reduction both urban rhizospheric soil as compared with urban peripheral soil (p = 0.020, Figure 5B) and in rural rhizospheric soil as compared with rural peripheral soil (p = 0.014, Figure 5C).The outcomes of a LEfSe indicated the potential prominence of Xanthomonadaceae within the rhizospheric zone (p = 0.007) (Figure 6).While functional differences were identified in bacterial communities in rhizospheric and peripheral soils, no statistically significant differences between the two soil types were identified using phenotype prediction across all of the following nine phenotype traits: aerobic, anaerobic, mobile element-containing, biofilm-forming, Gram-negative, Gram-positive, potentially pathogenic, and stresstolerant (Supplementary Figures S3, S4).
The results from the FUNGuild analysis, which categorizes fungal communities based on trophic mode, traits, and growth forms, reveal noticeable differences between the fungal communities inhabiting the rhizospheric and peripheral zones of the soil.Specifically, in peripheral soils, the trophic mode classification identified a comparatively heightened prevalence of pathotrophs (p = 0.042), coupled with an increased occurrence of plant pathogens as per the Guild classification (p = 0.011).Additionally, analysis of the growth form classification demonstrated significantly elevated proportions of smut (p = 0.026) and corticioid (p = 0.020) (Table 2).

Ascomycota fungal dominance in citrus orchard soils: ecological and agricultural implications
Ascomycota exhibited a substantial dominance within the citrus orchard soils.This pronounced presence of a specific taxon likely exerted a multifaceted impact on both the ecological landscape of citrus orchard soils and the agricultural productivity of the area.For example, the elevated prevalence of beneficial Trichoderma spp.(0.3%) would suppress plant pathogens.The principal mechanisms not only include mycoparasitism, antibiosis, and competition for resources and space but also involve the induction of resistance pathways in plants, leading to increased plant growth and nutrient uptake.Additionally, these fungi produce a diverse array of antifungal enzymes, such as chitinases and beta-1,3 glucanases.These enzymes exhibit synergistic effects with each other and with other antifungal enzymes and materials (Harman, 2006;Bae et al., 2011).Finally, rhizospherecompetent T. harzianum strain, such as T22, have been shown to promote root growth in a variety of plants (Harman, 2000).Conversely, the heightened abundance of potential pathogenic fungi, such as Fusarium spp.(0.4%), might increase the risk of disorders like root rot, consequently impeding plant growth and leading to diminished yield (Buttar et al., 2024;Kamble et al., 2024).
In terms of fruit quality, the prominent occurrence of Penicillium (0.4%) and Aspergillus spp.(0.4%).might lead to surface mold formation on citrus fruits, thereby influencing their visual appearance and gastronomic value.Moreover, certain strains within the Penicillium spp.could generate deleterious toxins, thereby posing latent risks to human health.Additionally, the augmented abundance of fungi similar to Chaetomium (0.4%) and Chaetomella spp.(0.3%), in turn, could substantively contribute to organic matter decomposition.
In the realm of farm management, a comprehensive understanding of the ecological roles of each species remains paramount.This comprehension is instrumental in devising judicious soil management strategies and disease control measures, designed to harness the proactive potential of beneficial species while mitigating the detrimental effects of harmful ones.

Deciphering urban-rural influences on citrus orchard soil: comprehensive analysis from diversity disparities to functional pathway differences
In the realm of diversity analysis, distinct differences were evident in the alpha diversity of bacterial communities between the citrus orchard soil of village A and those of villages B (p < 0.001), C (p < 0.001), and D (p < 0.001).The substantial disparity in species diversity became apparent through analyses of rarefaction curve plots, beta diversity plots, and taxonomic classifications at the phylum level (Figure 1).Considering the alignment of sampling design and realworld environmental observations, we inferred that the proximity to urban transportation zones played a pivotal role in driving differences in soil bacterial diversity between urban soil and rural soils.This inference was corroborated by the results of the functional disparity analysis of bacterial communities between urban soil and rural soil.
In the urban soil, the significant increase in the relative abundance of Formaldehyde assimilation II (RuMP Cycle) and formaldehyde oxidation I possibly stems from pollution sources, including vehicular emissions near urban areas, prompting bacterial communities to adapt metabolically to formaldehyde (Klein et al., 2022).Similarly, the notable increase in the relative abundance of toluene degradation V (aerobic) (via toluene-cis-diol) might be attributed to the prevalence volatile organic compounds, such as toluene, in urban environments, leading to metabolic adaptation among bacterial communities (Gelencsér et al., 1997).The significant rise in the relative abundance of anaerobic aromatic compound degradation (Thauera aromatica) is consistent with the presence of substantial aromatic compounds in the urban environment and the adaptation of bacterial communities to degrade these compounds (Keyte et al., 2016;Zhang et al., 2022).The urban orchards were closely situated alongside a motorized road, resulting in heightened vehicular activity during routine periods.Research on the impact of highway-related activities on soil microorganisms has shown that these activities exacerbate vegetation degradation along the road, significantly alter soil physicochemical properties, and cause heavy metal pollution, thereby affecting the diversity and community structure of soil bacteria.Furthermore, this disturbance exhibits a gradual increase in intensity with proximity to the highway (Liu Z. et al., 2021).In comparison of bacterial species in village A (urban) and the other villages (rural), Bacillus muralis, Bacillus fumarioli, and Bacillus flexus all showed significant increases in urban.Bacillus spp.possess heavy metal tolerance due to the presence of chromium reductase and nitroreductase enzymes, which are involved in the chromium (VI) reduction pathway (Luo et al., 2022;Shahraki et al., 2022;Ramli et al., 2023).This finding suggests that the changes in soil microbial community structure observed are influenced by traffic factors, consistent with the environmental factors and research results we have observed.It also suggests that in the study area, the application of B. muralis, B. fumarioli, or B. flexus for soil microbial remediation could be considered (Jebeli et al., 2017;Kulkova et al., 2023;Schommer et al., 2023;Yu et al., 2023).
Litter can increase the network complexity of soil microbial communities (Feng et al., 2022).Changes in litter and roots can affect soil microbial activities and nutrient cycling (Chen et al., 2022;Zhu et al., 2022).Potentially important physical differences between rural and urban areas were noticed during soil collection.We observed dispersed waste materials in rural orchards.These materials were primarily composed of discarded cigarette remnants, abandoned bags, garments, and packaging materials (Supplementary Figure S5), whereas similar materials were not observed in the urban orchards.Comparative analysis of abundant differential species between rhizospheric and peripheral soil.(A,B) Manhattan plot depicting the statistically significant differences in the abundance proportions of bacterial (A) and fungal (B) species between rhizospheric soil and peripheral soil.Patterns above the dashed line indicate sequences with statistical differences.(C) Heatmap displaying statistically significant bacterial species' differences between rhizospheric soil and peripheral soil, clustered by sample abundance.A p-value less than 0.05 was considered statistically significant.Dai et al. 10.3389/fmicb.2024.1403443Frontiers in Microbiology 10 frontiersin.orgBacterial distinctions in rhizospheric vs. peripheral soil revealed through LEfSe analysis.(A) The cladogram illustrates the microbial differences between rhizospheric soil and peripheral soil, highlighting only those microbial features that exhibit significant distinctions.Each circle represents a taxonomic level, such as phylum, class, and so forth.(B) The bar plot provides a detailed representation of the microbial differences between rhizospheric soil and peripheral soil.Correspondingly, there was a significant increase in the relative abundance of bacteria involved in nicotine degradation II (pyrrolidine pathway), indicating a potential bacterial metabolic adaptation to nicotine degradation.The significant elevation in the relative abundance of nylon-6 oligomer degradation might be linked to waste materials prevalent in rural environments.Fruit tree root system can affect soil microbial community structure, but it was not the main factor.According to our analysis, we observed a significant difference (p = 0.001) in the CMP-pseudaminate biosynthesis pathway between the rhizospheric soil and the peripheral soil (Figure 5).Pseudaminic acid, a sialic acidlike sugar, is involved in the synthesis of bacterial cell surface glycans, which in turn affects pathogenicity and adaptability through bacterial adhesion, invasion, and immune evasion mechanisms (Knirel et al., 2003;Hsu et al., 2006).The reduction in the CMP-pseudaminate biosynthesis pathway may represent a plant-microbe interaction, where plant roots release certain compounds that influence bacterial metabolism and biosynthetic pathways in the rhizospheric soil, contributing to a protective effect.However, it is important to note that our study also revealed a relatively high false discovery rate (FDR = 0.415, Supplementary Table S10), which may indicate a certain degree of false positives in multiple comparisons.Further research is needed to validate these findings and explore the biological significance of the observed differences in more depth.
These variations in bacterial functional pathways support our hypotheses regarding the impact of urban environmental factors on citrus orchard soil.The exogenous factors such as urban traffic and rural litter did have a significant impact on the microbiome in citrus orchards.These factors played a significant role in the substantial differences in bacterial diversity, potentially further influencing citrus tree growth and fruit yield.

Implementing practical strategies for citrus orchard soil management: insights and recommendations
By elucidating the microbial dynamics in citrus orchard soil during the crucial fruit maturation period, our findings can inform management strategies aimed at enhancing citrus yield and fruit quality.The identification of beneficial soil microorganisms capable of promoting citrus growth and disease resistance can lead to the development of targeted microbial-based biofertilizers and biocontrol strategies (Scott et al., 2006;Das et al., 2022;Wang et al., 2023), reducing reliance on chemical inputs and fostering long-term sustainability (Ali et al., 2019;Li et al., 2019;Zhou et al., 2020).
We propose the implementation of physical barriers in rural areas to reduce the quantity of pollutants such as cigarette butts and nylon.Urban orchards in city areas should be located far away from densely populated regions with high traffic density (Hagler et al., 2012;Ranasinghe et al., 2019).However, considering the current pollution situation, microbial measures can also be employed for improvement.For example, the Trichoderma harzianum T22 strain can be cultured on potato dextrose agar medium, and the spores can be collected for seed coating (Harman and Shoresh, 2007;Alinc et al., 2021) help improve plant growth and crop productivity under various abiotic stresses, including heavy metals and drought (Anbuganesan et al., 2024).In addition, synthesizing a biodegradable carbon nanoparticle from bacterial biofilm is also an option (Li et al., 2024;Savadiya et al., 2024).A more direct approach is to consider inserting genes encoding antifungal protein internal chitinase or external chitinase into the citrus genome alone or in combination (Bolar et al., 2000(Bolar et al., , 2001)).By utilizing amplicon sequencing as a non-invasive and costeffective method, we were able to gain insights into the microbial diversity and functional potential in citrus orchard soil at an unprecedented scale.Our study confirmed the beneficial application of metagenomics in addressing agricultural and environmental concerns.This included the identification of two pollutants that indeed affected microbial community functionality and proposed specific microbial methods for improving land heavy metal pollution and suppressing plant pathogens based on real-world microbial community structures.Consequently, future research could continue to utilize amplicon sequencing across different seasons and larger geographical scales, and shotgun sequencing methods could be employed to further validate the functionality of these microbes.
While our study has provided valuable insights into the microbial communities of citrus orchard soils, several limitations need to be considered.Firstly, our research solely relied on functional prediction methods to examine microbial community functional changes, which resulted in a lack of sufficient validation and interpretability.Future studies should incorporate more comprehensive functional analyses, such as employing shotgun sequencing methods to investigate specific differentially expressed genes.Moreover, our study did not include intervention studies using specific species, which could have strengthened the credibility of the proposed measures for soil health and agricultural productivity.Additionally, although we consider the fruit maturation stage as a crucial phase for investigation, it is essential to study soil in other growth periods to ensure the stability assessment of our findings.To maximize the comparability of sampling points, we applied specific restrictions, focusing on citrus orchards situated near rivers and lacking human management.Consequently, in urban settings, other suitable family farms meeting these criteria were not found.Further urban sampling and exploration across other soil environment categories are essential to enhance the universality of our conclusions.Finally, our study lacked a comprehensive assessment of environmental factors, which may have led to the omission of certain variables.Addressing these limitations in future research could further enhance our understanding of microbial communities in citrus orchard soils.

Conclusion
In this study, we aimed to elucidate the microbial community structures within the rhizospheric zone and peripheral soils of citrus orchards by employing environmental sampling and amplicon sequencing techniques.Our findings have shed light on the intricate interplay of urban environmental effects and root systems, both of which exert influence on the soil microbiota.Apart from the removal of specific pollutants, the application of B. muralis, B. fumarioli, and B. flexus to alleviate traffic pollution, and the use of Trichoderma spp.for plant pathogen suppression were considered viable solutions.The knowledge acquired from this study can be employed to optimize agricultural practices, augment citrus productivity, and foster sustainable agriculture in Zhejiang Province and beyond China.
Overall, our study, utilizing amplicon sequencing methods, identified that urban and rural environmental effects indeed influenced the microbial community structure in citrus orchard soils.A valuable approach for enhancing the stability and feasibility of conclusions in the future would involve the longterm observation of environmental and agricultural issues across various stages of plant growth using metagenomic methods.

FIGURE 1
FIGURE 1 Geographic distribution and taxonomy annotation of collection sites.(A) The specific locations of the four villages in the soil sampling process are marked in the diagram.The red arrows in the image represent the direction of the river's flow.(B,C) Stacked bar chart showing the average abundance of bacterial species (B) and fungal species (C) at the "phylum" level in the four villages.The chart highlights the top 8 bacterial species with the highest abundances.(A) Obtained from Google Satellite and annotated (https://www.google.com/maps/).

FIGURE 2
FIGURE 2 Comparison of species composition differences among different villages through alpha and beta diversity analysis.(A,B) Dilution curves of 16S (A) and ITS (B) sequences during the resampling process.(C,D) The alpha diversity differences of bacteria (C) and fungi (D) among the four villages were assessed using the Abundance-based Coverage Estimator (ACE) index.Each point in the box plot represents the ACE index of a sample.The lowercase letters "a" and "b" above the box plots indicate whether there is a statistical difference.The same letters indicate no statistical difference, while different letters indicate a statistical difference.(E) PCoA analysis based on Bray-Curtis distances between pairs of bacterial samples.The two dimensions with the highest explanatory power were plotted on the coordinate axes.(F) Stacked bar chart showing the average abundance of bacterial species at the class level in the four villages.The top 7 bacterial species with the highest abundances were highlighted.A p-value less than 0.05 was considered statistically significant.

FIGURE 3
FIGURE 3Phylogenetic trees of high-abundance species in rhizospheric and peripheral soil.(A,B) A maximum likelihood tree constructed from representative sequences of bacterial (A) and fungal (B) species with an average abundance proportion greater than 0.2% across all samples.The thickness of branches represents the magnitude of bootstrap values, and only values within the range of 0.2 to 1 are labeled.

FIGURE 5
FIGURE 5Differential bacterial species pathways in rhizospheric soil.(A-C) The bubble plot illustrates the differences in functional pathways of bacteria between rhizospheric soil vs. peripheral soil (A), urban rhizospheric soil vs. urban peripheral soil (B), rural rhizospheric soil vs. rural peripheral soil (C), presenting the top 10 pathways with the smallest p-values in ascending order.
for further analysis.

TABLE 1
Kruskal-Wallis test for differential ACE index comparisons between the bacterial population in soil from four villages.

TABLE 2
Selected results with statistically significant differences from guild classification comparisons of the microbial populations of rhizospheric and peripheral soils.
. Alternatively, soaking the seeds in a suspension containing Bacillus spp.for 2 h can